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Current evidence for dark matter in the universe does not exclude heavy composite nuclear-density objects 
consisting of bound quarks or antiquarks over a significant range of masses. Here we analyze one such pro- 
posed scenario, which hypothesizes antiquark nuggets with a range of B ~ 10 24 ~ 30 with specific predictions for 
spectral emissivity via interactions with normal matter. We find that, if these objects make up the majority of 
the dark matter density in the solar neighborhood, their radiation efficiency in solids is marginally constrained, 
due to limits from the total geothermal energy budget of the Earth. At allowed radiation efficiencies, the number 
density of such objects can be constrained to be well below dark matter densities by existing radio data over a 
mass range currently not restricted by other methods. 



I. INTRODUCTION 

Many different forms of astrophysical and cosmological ev- 
idence point to the existence of weakly interacting matter in 
an unknown form in the universe (see UJ for a recent re- 
view), with a total mass of order five times the inventory 
of what is observed as normal matter - gas, stars, and dust 
- giving a dark matter density estimated to be in the range 
0.4 — 1 GeV cm~ 3 in the solar neighborhood In addi- 

tion, the current best models for the early universe give strin- 
gent constraints on the content of interacting baryons in this 
dark matter at the time of nucleosynthesis JTJ. 

The rather compelling arguments for non-baryonic dark 
matter have led to a wide variety of efforts, both theoret- 
ical and experimental, to either postulate or directly detect 
new particles, beyond the standard model, that would satisfy 
the dark matter characteristics. However, there remain sev- 
eral "standard-model" candidates for dark matter, which, if 
not currently favored, have not yet been excluded over all the 
possible range of parameters. In particular, very massive ob- 
jects (at least compared to the <TeV scale of typical particle 
candidates for dark matter) can still satisfy the astrophysical 
constraints on dark matter if their masses are sufficient that 
the flux in typical detectors is extremely low, but not so large 
that they are excluded due to galaxy dynamics or gravitational 
lensing observations. This may be translated into a constraint 
on their interaction cross sections per unit mass; current limits 
require a/M < 0.1 cm 2 g _1 J5][6), a constraint that is in gen- 
eral easily satisfied by neutral objects with nuclear densities. 

Such objects must not interact with normal matter via 
strong or electromagnetic channels at the time of nucleosyn- 
thesis. One candidate is the so-called quark nugget, or 
strangelet, hypothesized originally by Witten 00, and de- 
veloped in many variations since then. Quark nuggets can be 
neutral and metastable at their formation during the quantum 
chromodynamics phase transition of early-universe evolution, 
and thus do not undergo significant further interactions at nu- 
cleosynthesis, therefore evading the constraints on baryonic 
content. 

Although still a matter for debate, the possibility of quark 
nugget color superconductivity (9), in which quarks near the 
Fermi surface of the nugget form correlated Cooper pairs, fa- 
vors their possible stability ifTol . One of the more attractive 
aspects of these objects as candidates for dark matter is that 



the physics of their formation and interactions is in principle 
calculable according to the standard model, although such cal- 
culations can in practice be prohibitively difficult. 



II. ANTIQUARK NUGGETS 

A recent novel application of this model lfl6lTT9ll postu- 
lates that both antiquark and quark nuggets are formed in the 
early universe with a ratio of 3:2, and the current observed 
baryon matter-antimatter asymmetry arises only because the 
antibaryons are hidden in the excess of antiquark nuggets 
(AQN), which, along with the quark nuggets form the bulk of 
the dark matter. AQN have the same kinetic energy as normal 
quark nuggets, and the transfer of this energy may be observed 
in seismic [ 12] or thermal events produced in the Earth's crust. 
However, in addition to kinetic energy transfer, AQN sweep 
up and annihilate with normal matter along their track, lead- 
ing to potentially much more energetic signatures and much 
higher rates of radiative energy deposition. 

Current constraints from seismic energy deposition in the 
Earth and Moon indicate that quark nuggets can only satisfy 
dark matter density for baryon numbers B < 7 x 10 28 [ 14 1 (but 
see the possible detection of an event in ref. |[T3"1 . and the 
response in [15]). Limits from non-detection of compatible 
events in the Lake Baikal detector J20l require B > 1 .2 x 10 20 . 

Taking an approximate geometric mean value of these con- 
straints, a baryon number of B > 10 24 (~ 1 .6 gm) gives a flux 
of AQN at Earth, assuming they are virialized with Galactic 
velocities of order 200 km s _1 , of order several per km 2 per 
year if all AQNs were close to this mass; actual fluxes will 
depend of course on the assumed mass spectrum in the solar 
neighborhood. 



A. AQN thermal emission 

It is instructive to consider the flux of AQN at this 
mass scale to determine the rate of energy deposited in the 
Earth. Recent detailed calculations of the emissivity of AQN 
when accreting normal matter have been carried out using a 
Thomas-Fermi model lfT6llT7ll . In this case the spectral emis- 
sivity of a nugget at effective temperature T for photon ener- 
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lies well below the electron mass m e is given by: 
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where Osb is the Stefan-Boltzmann constant, a is the fine- 
structure constant, and CO is the angular frequency of the radi- 
ation, and 
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(Here ratios such as co/T are given with implicit values of 
Planck's and Boltzmann's constants h, k where necessary to 
rationalize the units.) This result, converted to spectral power 
dP/d(0 — AizR^dF / d(i>, is plotted in Fig.[T]for the case of an 
AQN of B = 10 24 , and a nugget radius estimated flSE) to be 
R„ = 10 m, for three different AQN surface temperatures. 

The integral emissivity over all frequencies, determined 
from this spectral emissivity is 
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and is typically a fraction of order 10~ 6 of blackbody emis- 
sion ITBirnn . This power is however distributed with a spec- 
tral emissivity very different from a blackbody, since it is 
nearly flat at low frequencies as equation[JJand Fig.[JJshow. 
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FIG. 1: Spectral density of emission from a B ~ 10 24 antiquark 
nugget for three values of the effective surface temperature. 



nugget's path, by normalizing the total emission (the integral 
of this equation over all co) to equal the fraction 1 — g of an- 
nihilation energy that the nugget thermalizes from its accre- 
tion IfTBI . thus F tot = (1 — g)F ann . For atmospheric and sur- 
face densities encountered at Earth, the density-dependent ef- 
fective temperature is expected to be IfTBI 
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where / is related to the accretion efficiency of the AQN, u n 
is the AQN velocity relative to Earth, and p is the density of 
the medium. The annihilation efficiency terms are somewhat 
uncertain; typical values given are g ~ 0. 1 and / ~ 0.067 IfTBI . 

A firm upper bound on the maximum AQN temperature 
Tmax comes from the requirement that the radiation pressure 
just outside the surface of the AQN should not exceed the 
level required to accelerate the accreting matter away from 
the path of the nugget. A previous qualitative analysis |fT9l 
of the effects of radiation pressure of incoming material in 
the Earth's atmosphere estimated that T would saturate at 
densities p ~ 10~ 3 g cm~ 3 in the lower atmosphere, giv- 
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10 keV. This estimate did not account for the 



likely photo-ionization of the region very near the nugget and 
thus probably underestimated T majc ETI I. since the radiation 
pressure primarily depends on the cross section for photon 
momentum transfer with incoming nuclei. Another possible 
accretion-limiting effect was also discussed in reference |fT9l , 
that of momentum transfer to the incoming nuclei via scat- 
tering off positrons in the AQN electrosphere surrounding the 
nugget. We first discuss the radiation pressure bound, the re- 
turn later to the effects of positron scattering. 

To illustrate these effects, consider the case of a nugget with 
B = 10 24 , and T ~ 10 keV (~ 10 8 K) accreting material at 
sea-level Earth-atmospheric density p ~ 10~ 3 g cm~ 3 at a ve- 
locity of 200 km s~'. For this case, the AQN luminosity is 
Ltot = 4%RlF tut = 2 x 10 6 W. The intensity just outside the 
surface of the nugget is I sur f = F tot ~ 1.6x 10 19 W m 2 . Since 
a large fraction of this intensity is emitted as soft X-rays, for 
which the photo-ionization cross section on air will far ex- 
ceed the momentum-transfer cross section On, the incoming 
material will be fully ionized before it approaches the AQN 
surface. The residual momentum-transfer cross section on the 
stripped nuclei is not well-documented, but it cannot exceed 
the incoherent atomic scattering cross section at keV photon 
energies, thus we estimate On ^5 10~ 2 barn = 10~ 30 m 2 . The 
required acceleration to displace the nuclei from the oncom- 
ing nuggets path is of order a m ,„ = Am/A? ~ u%/R„, and thus 
the radiation pressure is determined from 
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B. Temperature upper bound 

The temperature of the AQN in equation [JJ is determined 
from the amount of matter accreted and annihilated along the 



where mjy ~ 2.2 x 10 26 kg is the typical atomic mass, for 
nitrogen in this case. This yields, for the example above, 
I surf < 3 X 10 30 W m 2 , many orders of magnitude above the 
the AGN radiance in this case. Bounding the surface emissiv- 
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ity with this then bounds the AQN temperature: 
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and since the AQN radius depends on the number of baryons 
BasR„ - KT^fl/lO 24 ) 1 / 3 m 2 H2): 



r mra < 4.5 MeV 
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Of course, at several MeV AQN surface temperatures, exter- 
nal nuclear interactions and associated energy release will be- 
come important, and the luminosity of the nugget will be so 
high that other limiting effects are likely to obtain well before 
these temperatures are reached. 

Another potential limiting effect is that the ram pressure 
of the fraction of stripped nuclei that mechanically collide 
with positrons in the electrosphere of the nugget will form 
a bow shock to the propagating nugget, and this will limit the 
ingress of other matter. This constraint was discussed briefly 
by Lawson |fT9ll , where it was argued that once the kinetic 
temperature of the positrons approached that of the incom- 
ing nuclei, they would begin to scatter off the positrons and 
a negative feedback condition would obtain. However, the 
momentum transfer of the electrosphere to incoming nuclei is 
given by the Rutherford scattering cross section for positrons 
on the stripped nuclei. As the nugget temperature increases, 
the Rutherford scattering cross section decreases quadratically 
with the effective temperature Or(T) °c T~ 2 . In contrast, the 
AQN electrosphere temperature increases only slowly with 

4/17 . . 

annihilation rate T 

oc F ann 111611 , causing this process to de- 
crease in efficiency in deflecting nuclei as the accretion pro- 
cess increases. Thus it appears that positron scattering cannot 
produce a negative feedback accretion-limiting condition. 

A more likely source of negative feedback may come from 
backscattering at the quark-matter interface of the nugget; the 
accretion efficiency term / used above arises from this pro- 
cess iTTBI . but it is unclear what the temperature dependence 
of this effect might be. An estimate of this process is beyond 
our scope, so in what follows we assume that AQN surface 
temperatures approaching 100 keV (a few percent of the ra- 
diation pressure T max derived here) in solid materials are not 
excluded yet by any accretion constraints. 



III. AQN INTERACTING AT EARTH 

A. Lithosphere interactions 

We first consider interactions of these AQN in Earth's litho- 
sphere. The kinetic energy loss for either quark or anitquark 
nuggets is given by lfl2l 
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where A„ is the cross sectional area of the nugget, and prx) 
is the density along the track x. Using the average density 

pW = (p), 
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If an antiquark nugget stops in the Earth through loss of ki- 
netic energy, it will then subsequently completely annihilate, 
and the energy deposition in that case will equal its remaining 
mass energy. For AQN to remain viable dark matter candi- 
dates, the total power contribution through this process to the 
Earth's thermal energy budget must not exceed that of known 
sources. The current geothermal energy budget of the Earth 
is P geo ~ 44 ± 3 TW GU G2, and of order half of this must 
be attributable to radionuclide decay; the remainder is still a 
subject for debate, although a major fraction must be resid- 
ual heat from the gravitational collapse of formation. If we 
allow that of order P ge o/^ ~ H TW of the current geother- 
mal energy budget could be available to external heating from 
AQN annihilation, then the rate of captured AQN must satisfy 
dm/dt < P geo /(4c 2 ) = 0.12 grams/sec. At the current firm 
lower bound for AQN baryon number Z?„„„ ~ 10 20 =0.16 mg, 
a flux of AQN at this mass, equal to the dark matter density, is 
of order 10 6 s _1 over the whole Earth, thus the capture prob- 
ability must not exceed about 10~ 3 per nugget. The required 
flux to match the DM falls as B^ 1 , and the mass-energy rises 
as B, so this constraint is constant with AQN mass. So far 
we have ignored the energy deposited during the transit of the 
nugget; we return to that below. 

If we require that no more than 0.1% of all AQN lose 
enough kinetic energy via equation[7]to be captured by Earth, 
this translates into a requirement that the velocity attenuation 
in equation [8] above can only fall below the escape velocity 
in 0.1% of all AQN tracks. Taking m„(0) = 200 km s" 1 , 
and the Earth escape velocity u esc ~ 11 km s , we evolve 
an initial Maxwell-Boltzmann velocity distribution by assum- 
ing equation [8] above, with a mean travel distance of x — 
4Re/3, the mean chord distance through the Earth for ran- 
dom tracks 11221 1231 . By then requiring that the cumulative 
evolved Maxwellian have no more than 0.1% of its final ve- 
locities below u esc , we find the following constraint if all of 
the dark matter consists of AQN of mass equal to or greater 
thanB„„„: 
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where we have used a mean density of (p) = 5500 kg m 3 
for the Earth, and a nugget area A„ — %R 2 , with R„ = 
lO-^B/lO 24 ) 1 / 3 m. It thus appears that geothermal consider- 
ations rule out AQN of masses less than about 4 grams, well 
above prior constraints. This limit can only be evaded if the 
velocity distributions of the AQN are decidedly non-virial, but 
a similar constraint will obtain on whatever velocity distribu- 
tion is present. 

Now consider the emissivity in lithosphere transit of a flux 
of AQN well above B m j„ given here, with a nearly mono-mass 
spectrum with B ~ 10 25 . As the nugget enters the solid crust 
at 200 km s _1 , the temperature rises to around 120 keV (us- 
ing expected values for /, g above, and a mean density of 
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5.5 gm cm" 3 , giving an initial luminosity Pi = 471 R 2 F tot ~ 
3.7 x 10 11 W. At this mass, the velocity attenuation length 
is M/(A„(p)) ~ kRe, and using the mean chord of 4Re/3, 
the velocity at exit is Uf ~ (2/3)«„(0), so the mean velocity 
is of order 0.8m„(0), and thus the average power is 0.8P, ~ 
3 x 10 W. The integrated rate of energy being continuously 
deposited, if the dark matter consists entirely of nuggets of 
this mass, is P tot ~ 48 TW. This level of thermal energy also 
exceeds the current ~ 44 TW geothermal energy budget of the 
Earth ll24l |251 . Using again the requirement that AQN con- 
tribute no more than 1/4 of the current geothermal energy, 
we may place a constraint on the maximum temperature at 
B ~ 10 25 : in this case T max < 120 keV(ll/48) 4 / 17 = 85 keV. 

It might appear that the T 4 dependence of the AQN lu- 
minosity does not leave much headroom for larger nugget 
masses, since the accretion rate grows with the cross section 
of the nugget. Clearly, the radiation pressure constraint used 
in deriving equation|6]above is far less restrictive than the con- 
straints from geothermal power. However, the temperature is 
to first order only dependent on the density of the medium, 
rather than nugget mass. Also, since the accretion cross sec- 
tion grows as B 2 ' 3 , and the flux required to match the dark 
matter density decreases as B , higher masses are possible, 
but again will be marginally close to violation of the geother- 
mal constraints, unless the maximum temperature or accretion 
efficiency is lower than initial estimates. 



B. AQN meteors? 

Propagating in the Earth's lower atmosphere at a temper- 
ature of 10 keV, an AQN with B ~ 10 24 produces megawatt 
total bolometric luminosity. However, its apparent visual sig- 
nature in the optical band would not necessarily be dramati- 
cally different than a very fast meteor. At u„ ~ 200 km s" 1 it 
is only a factor of 3 faster than the fastest meteors, although 
its trail could be much longer. At this temperature, the AQN 
power in the optical band is of order 1 kW, equivalent to rela- 
tively bright meteor with visual magnitude of my ~ — 1 l26l . 
However, the nugget would not achieve this temperature un- 
til near ground level; at typical ~ 80 km meteor altitudes, the 
AQN temperature would be a factor of 10 lower, and the lu- 
minosity more like 100 mW, thus making it invisible to the 
naked-eye at higher altitudes, where it would be much more 
widely visible. It is not hard to see why such events could 
evade casual observation. 



C. Thermal radio signatures 



the atmosphere is thus of order |[T9l 



While the visual emissivity of an AQN with B ~ 10 pass- 
ing through the terrestrial atmosphere may be easily missed, 
the flat spectrum displayed in Fig. [T] produces a surprisingly 
strong broadband radio signature, of order 10 mW GHz -1 in 
the VHF to microwave band even at the lowest AQN tempera- 
ture considered, and far more at the higher ones. The expected 
radio flux density in this spectral region for an AQN transiting 
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where we have ignored the weak logarithmic dependence. For 
an AQN transiting a solid material, viewed from outside the 
material, there are additional terms due to the attenuation of 
radio signals in the solid material, and the Fresnel transmis- 
sion coefficient for the emission as it passes through the 
interface. Thus 
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where d is the pathlength of the radio emission in the solid, 
and Latten is the field attenuation length in the medium. 

Receiver thermal noise levels at a typical receiver system 
temperature of T sys < 300K, by comparison, are typically 
P n = kT sys Af < 5 pW in a GHz of bandwidth. The broad- 
band AQN radio flux density is thus likely to be well above 
thermal noise for large distances from the track. This will of 
course depend on the mean distance (D) from the track, as 
well as the time x over which the AQN track remains in the 
primary field-of-view, or half-power beamwidth H, of a given 
antenna. This in turn depends on the antenna gain (or direc- 
tivity) G. For moderately directive antennas of a few dBi of 
gain or more, G ~ 27000 (//°)" 2 where the beamwidth H is 
given in degrees here, and thus for u n ~ 200 km s , G ~ 10, 
and (D) ~ 100 km, the in-beam residence time is 0.3-0.5 sec- 
onds; however, receiver gain instabilities make it practical to 
limit the integration time X ^ 0.1 seconds, giving several sam- 
ples per transit per beam. An antenna of constant gain G over 
a passband from fa to fa has an average effective area (for a 
flat spectrum source) of A e ff = 2Gc 2 / [Anfafa) and the min- 
imum detectable signal power of this antenna with receiver 
bandwidth A/ = fa — fa and integration time X is 
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Assuming the integration time x is matched to the expected 
beam crossing time for u„ ~ 200 km s -1 , the limiting sensi- 
tivity is 
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Comparing this to equations 10 and 11 indicates that such 



events are detectable with a modest antenna collecting area 
and receiver out to distances of several hundred km, even at 
the lowest AQN temperature considered here. The advantage 
in detection of thermal radio emission as compared to other 
possible forms of beamed emission, such as geosynchrotron 
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emission considered in reference lfl9l . is that the acceptance 
solid angle does not require observation very close to the axis 
of the AQN velocity. Thus for an isotropic flux of AQN, ther- 
mal radio detection will have a far greater probability. 
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FIG. 2: Estimated flux density and sensitivity for ANITA for two dif- 
ferent AQN temperatures and three different distances, for AQN tran- 
siting in atmosphere (T = 10 keV), or Antarctic ice (T = 80 keV), 
with typical estimated attenuation losses and Fresnel coefficient. The 
thermal noise levels for two different integration times are shown. 



IV. THERMAL RADIO DETECTION PROSPECTS 

Given that the distance scale for detection is plausibly hun- 
dreds of kilometers even at the lowest AQN temperatures, 
and perhaps much more at higher temperatures, it is evident 
that a ground-based detector is at a disadvantage compared to 
the synoptic field-of-regard of suborbital or orbital platforms. 
Since the AQN luminosity appears to rise by orders of mag- 
nitude once it enters solid material, and since terrestrial ice in 
the Earth's cryosphere is highly transparent in the VHF and 
UHF radio range ll27l l28l . suborbital/orbital observations of 
Antarctic or Arctic ice sheets afford perhaps the most sensi- 
tive possible channel for AQN detection. We thus conclude 
by estimating to first-order the sensitivity of this approach, 
using the parameters of the Antarctic Impulsive Transient An- 
tenna (ANITA) suborbital pay load |29|, which has completed 
two flights (ANITA- 1: 2006-2007; ANITA-2: 2009-2010) 
and is scheduled to complete a third next year. ANITA has 
enough directional capability to establish the velocity of an 
AQN candidate, a crucial discriminator against other possible 
background events such as meteors. 

For a radio detector viewing Antarctica synoptically from 
stratospheric altitudes, as in the case of ANITA, the horizon 



is at a distance of 600 — 700 km, and the area viewed is over 
1M km 2 out to the horizon. To illustrate the range of sen- 
sitivity, Fig. [2] plots the AQN signal and thermal noise lev- 
els based on equations 10 11 & 13 above, using instrument 
parameters for ANITA-2 and ANITA-3 [30|, for a range of 
AQN masses that are currently unconstrained, for several dif- 
ferent distances. For ANITA-2, the x = 1 /js integration time 
arose in the ambient radio-frequency (RF) power monitoring 
system, an auxiliary detector to the primary 2.6 Gsample/sec 
waveform recorder which captures only a 100 ns time window. 
The ambient RF power monitor samples each antenna signal 
at about 8 Hz, but with a much shorter integration time due 
to analog-to-digital conversion related constraints. However, 
these samples occur for both polarizations, and there are 2-3 
antennas sampling each azimuthal direction at 22.5° azimuth 
intervals. Since the thermal AQN signals are unpolarized, the 
two ANITA polarization samples are independent, as are the 
2-3 different antennas per azimuthal sector, and thus an AQN 
signature can be detected with high-confidence by requiring 
the signals to appear in a majority coincidence of all of these 
independent channels. For ANITA-3 the design will allow for 
much longer integration times per RF power monitoring sig- 
nal, and thus the sensitivity will be substantially improved. 

It is evident from Fig.|2]that AQNs transiting either the at- 
mosphere or ice sheets can be detected, but to ensure that these 
are distinguished from the many forms of anthropogenic and 
other radio-frequency interference, the AQN track needs to be 
detected over an azimuth range such that a velocity can be un- 
ambiguously determined. This requires a projected azimuthal 
span for the track of order A(j) > 20° for ANITA (29), so that 
at least two adjacent azimuthal sectors show a signal; this is 
possible over this range of azimuth since the ANITA antenna 
beams overlap each other in adjacent azimuthal sectors. 

To estimate the limiting sensitivity for ANITA, we must in- 
tegrate over the acceptance area, solid angle, and detection 
efficiency. The number of detected events N as a function of 
the baryon number B of the AQN can be written 
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where T b s is the integrated observation time, is 
the distance to the horizon, 7 ifi) is the flux of AQN, 
and T,(B,r, 8,(j),p) is the instrumental detection efficiency 
(bounded between [0,1]) as a function of baryon number, dis- 
tance r, angular track directions 0,(|), and medium density p. 

The requirement on azimuth span above constrains the 
track to be relatively horizontal, such that at distance D 
from the payload, the zenith angle range is limited by 
A8 ~ (10 km) / (Di/«A(j)) for events in the atmosphere, A9 ~ 
(2.0 km) / (DsinAty) for in-ice tracks. Using these constraints, 
and requiring that the signal- to-ratios Saqn^mos/^s > 3 or 
SAQN,ke/<5s > 3 for a trigger, we have simulated the detec- 
tion of AQNs with ANITA to determine the limiting sensi- 
tivity using Monte Carlo methods. The results are shown in 
Fig. [3] where we have shown the curves based on the flown 
ANITA-2, and planned ANITA-3 instruments. For ANITA- 
2, this applies to data that has been already acquired, but has 
not yet been analyzed for this type of signal; here we assume 
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FIG. 3: Existing and projected limits on quark nugget fluxes with 
respect to dark matter density ( diagonal hatched curve ). The ANITA 
curves are estimates from this work. 

that no background events are found. For ANITA-3 we use 
projected performance estimates provided by the ANITA col- 
laboration, and a 30 day livetime is assumed, similar to what 



was achieved for ANITA-2. The inflection in the sensitivity 
curve in each case is due to the crossover of the effects of 
ice detection, which is more efficient at lower AQN masses, 
and atmospheric detection, which has a larger available solid 
angle and is more effective at higher AQN masses. For both 
ANITA-2 and ANITA-3, the sensitivity eventually saturates 
the available area for very high AQN masses, and the curves 
flatten out. 

We have included in Fig. [3] the constraint from geothermal 
power derived above, and we note the other two actual (rather 
than projected) limits are from the Lake Baikal detector 11201 , 
and from analysis of the Lunar seismic noise detected in the 
Apollo 11 mission lfl4l . For the IceCube curves, both are 
still projections, although the data for IceCube 22 already ex- 
ists. For IceCube 80, the projection is for three years of data, 
which will be achieved in early 2014 1311 . Analysis of exist- 
ing ANITA-2 data for these event signatures has begun, and 
results may be expected within the next year. Thus it appears 
that this very interesting region of parameter space for quark 
nugget dark matter is within reach of several current experi- 
ments, and we can hope for either detections or compelling 
limits in the near future. 
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